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Abstract 

q 

It has been recently shown that the quantum Boltzmann equations may be relevant for 
the leptogenesis scenario. In particular, they lead to a time-dependent CP asymmetry 
which depends upon the previous dynamics of the system. This memory effect in the 
CP asymmetry is particularly important in resonant leptogenesis where the asymmetry 
is generated by the decays of nearly mass-degenerate right-handed neutrinos. We 
study the impact of the non-trivial time evolution of the CP asymmetry in the so- 
called Minimal Lepton Flavour Violation framework where the charged-lepton and the 
neutrino Yukawa couplings are the only irreducible sources of lepton- flavour symmetry 
breaking and resonant leptogenesis is achieved. We show that significant quantitative 
differences arise with respect to the case in which the time dependence of the CP 
asymmetry is neglected. 



X 



1 Introduction 



The observed baryon number asymmetry of the Universe t\b = ns/n-y = (6.3 ± 0.3) x 1CP 10 [T] can 
be explained by the mechanism of thermal leptogenesis [2] [3j 0], the simplest implementation of 
this mechanism being realised by adding to the Standard Model (SM) three heavy right-handed 
(RH) neutrinos, i.e. in the framework of the seesaw (5]. In thermal leptogenesis the heavy 
RH neutrinos are produced by thermal scatterings after inflation and subsequently decay out- 
of-equilibrium in a lepton number and CP-violating way, thus satisfying Sakharov's constraints 
(HI [7] . The dynamically generated lepton asymmetry is then converted into a baryon asymmetry 
due to {B + L)-violating sphaleron interactions |8j. 

If RH neutrinos are hierarchical in mass, successful leptogenesis requires that the RH neutrinos 
are heavier than about 10 9 GeV |9j. Hence, the required reheating temperature after inflation 
cannot be much lower [3] [10] [11] [12] . In supersymmetric scenarios this may be in conflict with the 
upper bound on the reheating temperature necessary to avoid the overproduction of gravitinos 
during reheating |13j . In the resonant leptogenesis scenario [14 this tension may be avoided. If 
the RH neutrinos are nearly degenerate in mass, self-energy contributions to the CP asymmetries 
may be resonantly enhanced, thus making thermal leptogenesis viable at temperatures as low as 
the TeV. Resonant leptogenesis seems therefore a natural possibility. 

Nearly degenerate RH neutrinos naturally arise in the context of models satisfying the hy- 
pothesis of Minimal Flavour Violation (MFV) |15[ [T6] [T7] . In the quark sector, where the MFV 
hypothesis has been formulated first, the MFV ansatz states that the two quark Yukawa couplings 
are the only irreducible breaking sources of the flavour-symmetry group defined by the gauge- 
kinetic Lagrangian [16 . In generic models satisfying this hypothesis, quark Flavour Changing 
Neutral Currents (FCNC) are naturally suppressed to a level comparable to experiments and 
new degrees of freedom can naturally appear in the TeV range. The extension of the MFV 
hypothesis to the lepton sector (MLFV) has been formulated in Ref. [T7] (and further studied 
in Refs. |18[ [T9] I20|). where a number of possible scenarios, depending on the field content of 
the theory, have been identified. The case more similar to the quark sector and more inter- 
esting from the point of view of leptogenesis is the so-called extended field content scenario of 
Ref. |17) . Here the lepton field content is extended by three heavy RH neutrinos with degenerate 
masses at the tree level. The largest lepton-flavour symmetry group of the gauge-kinetic term is 
G max = SU(3)e x SU(3) e x 0(3) at and, according to the MLFV hypothesis, it is assumed that this 
group is broken only by the charged- lepton and neutrino Yukawa couplings A e and X u . In relation 
to leptogenesis, the key feature of this scenario is that the degeneracy of the RH neutrinos is lifted 
only by corrections induced by the Yukawa couplings, so that we end up with a highly constrained 
version of resonant leptogenesis. Within this set up, the viability of leptogenesis has first been 
considered in the one-flavour approximation [21] and subsequently in the flavoured case |22| [23] . 

Resonant leptogenesis in models satisfying the MLFV hypothesis is the subject of the present 
paper. Our analysis turns out to be necessary in view of the recent results achieved by studying 
the dynamics of thermal leptogenesis by means of quantum Boltzmann equations |24| (for an 
earlier study, see Ref. [25]). Quantum Boltzmann equations were obtained starting from the 
non-equilibrium quantum field theory based on the Closed Time-Path formulation and have an 
obvious interpretation in terms of gain and loss processes. However, they differ from the classical 
Boltzmann equations since there appear integrals over the time where theta functions ensure that 
the dynamics is causal. In the classical kinetic theory the scattering terms do not include any 
integral over the past history of the system which is equivalent to assume that any collision in the 
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plasma does not depend upon the previous ones. On the contrary, quantum distributions possess 
strong memory effects and the thermalization time obtained from the quantum transport theory 
may be substantially longer than the one obtained from the classical kinetic theory. Furthermore, 
and more importantly, the CP asymmetry turns out to be a function of time, its value at a 
given instant depending upon the previous history of the system. If the time variation of the 
CP asymmetry is shorter than the relaxation time of the particles abundances, the solutions 
to the quantum and the classical Boltzmann equations are expected to differ only by terms of 
the order of the ratio of the time-scale of the CP asymmetry to the relaxation time-scale of the 
distribution. In thermal leptogenesis with hierarchical RH neutrinos this is typically the case. 
However, in the resonant leptogenesis scenario, RH neutrinos are almost degenerate in mass and 
the CP asymmetry from the decay of the i-th RH neutrino iVj is resonantly enhanced by the j-th 
neutrino if the mass difference (Mj — Mi) is of the order of the decay rate of the RH neutrinos. 
The typical time-scale to build up coherently the time-dependent CP asymmetry is of the order 
of (Mj — Mj) -1 [24, 26j, which can be larger than the time-scale for the change of the abundance 
of the RH neutrinos. This tells us that in the case of resonant leptogenesis significant differences 
are expected between the classical and the quantum approach. 

In this paper we perform a detailed study of the role of quantum memory effects in the resonant 
leptogenesis scenario within the MLFV hypothesis showing that memory effects may substantially 
change the prediction for the baryon asymmetry. We consider both the normal hierarchical (NH) 
and the inverse hierarchical (IH) cases for light neutrinos and also consider the role played by 
flavour in leptogenesis [27] 1281 121J1 130] . 

Our analysis is organized as follows. In Section 2 we provide a brief summary of thermal 
leptogenesis, focussing on the impact of the quantum memory effects. Section 3 contains the key 
ingredients of the MLFV scenario. In Section 4 we set the stage for the study of the relevance of 
non-equilibrium quantum effects in MLFV-leptogenesis. We then proceed discussing the impact 
of quantum effects on the two scenarios in which CP violation arises from the RH sector in Section 
5, and only from the PMNS matrix in Section 6. We give our conclusions in Section [7] 



2 A brief summary of thermal leptogenesis 

In order to set the stage for the subsequent analysis about the impact of non-equilibrium quantum 
effects in the MLFV framework, here we first briefly recall the general features of the thermal 
leptogenesis scenario, and later discuss the restrictions imposed by the hypothesis of Minimal 
Flavour Violation. 

We consider a model where three right-handed (RH) neutrinos N% (i = 1, 2, 3) with Majorana 
masses M3 > M2 > M\ and Yukawa couplings X u are added to the Standard Model (SM) 
Lagrangian. Working in the basis in which the Yukawa couplings for the charged leptons A e are 
diagonal, the Lagrangian density is given by 

C = C SM +(^M l Nf + N i (X u ) ia e a H + e a (X e ) a £ a H c + h. C }j , (1) 

where l a and e a indicate the lepton doublet and singlet with flavour (a = e, jx, r) respectively, 
and H is the Higgs doublet whose neutral component has a vacuum expectation value equal to 
v = 174 GeV. For heavy RH neutrinos (Mj ^> v), light neutrino masses are obtained via the 
see-saw (type I) mechanism |5j 

m u = U*m u U j = \lM~ 1 X u v 2 , (2) 
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where U is the PMNS mixing matrix (and we adopt the convention of placing a hat over di- 
agonal matrices with real and positive elements). In this framework the baryon asymmetry is 
generated by weak sphaleron processes converting the non-zero lepton number produced by out- 
of-equilibrium decays of the heavy RH neutrinos. 

The key quantities controlling the production of a net lepton number are the CP violating 
asymmetries in the iV, decay rates 



= F(Nj -> i a H) - F(Nj -> l a H) 
€ia ~ T(Ni ^ IJT) + TjNi - I a H) ' 

The inclusion of quantum effects discussed in the introduction (for technical details see Refs. E 
[26] ) introduces a time dependence in the CP asymmetry 
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where Tj = J2/3^(^j ~ * = (A^At)jjMj/(87r) is the total decay rate of the j-th RH neutrino, 
H is the Hubble expansion rate, z = Mx/T, T denotes the temperature, g s and g v are the self- 
energy and the vertex correction functions respectively and X{ = (Mi/M\) 2 (a short summary of 
the Boltzmann equations can be found in the Appendix A) . The combination of Yukawa couplings 
appearing in Eq. Q is quite constrained under the hypothesis of MLFV and, as we discuss below, 
the requirement of non- vanishing £j a leads to non-trivial constraints on the free parameters of this 
framework. 

In the quantum Boltzmann approach, the typical time-scale for the variation of the CP asym- 
metry is 
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1 2H(T) 2i?(M x ) ~ Mj-Mi AM i; ' ^ 

As a consequence, quantum effects are expected to be sizable if 1/AMji is larger than the time- 
scale for changing the abundance, 1/Tj. In other words, the oscillation frequency AMji has to 
be sufficiently smaller than Tj, so that the factors m do not effectively average to one. 
Under these conditions, the amplitude of the "sin" term in m^'^z) is also enhanced, which will 
turn out to be a crucial effect. If 1/AMji is smaller than the time-scale for changing the particle 
abundances, then the CP asymmetry may be averaged over many scatterings [M] and it reduces 
to the classical value ej Q = X^y* e a'^- 

The above discussion leads us to formulate a quantitative criterion for the importance of 
quantum effects, namely AMji ~ Tj- This criterion can be naturally satisfied if RH neutrinos are 
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nearly degenerate, as first pointed out in and therefore in models based on MLFV. So the next 
task is to identify the parameters controlling the ratio AMjj/Tj within the MLFV framework. 



3 The MLFV scenario 

The MLFV approach with extended field content [T7] is based on the assumption that the largest 
possible flavour symmetry group of the lepton sector, SU(3)i x SU(3) e x 0(3)n, is broken only 
by two irreducible sources: the Yukawa couplings Ag and A° transforming as (3, 3, 1) and (3, 1,3) 
respectively. 

In the limit of vanishing Yukawa couplings, the 0(3)n symmetry is exact and RH neutrinos 
are degenerate at the common scale M®. Switching on the Yukawa interactions A^ and Ag, the 
mass degeneracy is removed as an effect of those combinations of Yukawa couplings transforming 
as (1,1,6.) M - 



M v = M, 



1 + c« {h% + hf) + c? {{hlf + {hf?) 
+ 4 2) hlhf + cf hfhl + cf {h° e + hf) + .... 1 , (6) 



where h° v = A°A° f , h° e = A°A° f A°A° f and one can choos^j] Ag = rh e /v. The coefficients c's are 
arbitrary parameters smaller than 0(1). If one interprets them as arising from radiative effects, 
it is natural to have ~ l/(4-7r) 2 and cf^ ~ (c^) 2 . From the unitary matrix U diagonalising 
M u , one obtains M v and A^ of Eqs.O and rt2b: 



M U = UM U U T , X U = UX° U . (7) 

Notice that, due to the smallness of the c's, it is natural to expect U to be close to the unity 
matrix (or possibly a permutation matrix, given the ordering convention M\ < M2 < M3). 

One can see A^ as the neutrino Yukawa coupling associated to degenerate RH neutrinos. 
Accordingly, a convenient parameterisation to be exploited in the following is [31] 

\l= l - ^W v H ^ u U oj , (8) 





where m° = diag(m?, m°, m§), U° = R 2 3(e 23 )T s oR 13 (9%)Tl R 12 (e 12 )dmg(e i ^ , , 1), T s o 
diag(l,l,e i5 ) and H is an orthogonal hermitian matrix 



H = e iq> , $ = ( -0i 03 I • (9) 

■02 

Clearly, the smaller the splittings among the RH neutrinos, the closer rh v and U are to and 
U° respectively. Since the approximations rn v ~ m|J and U ~ U° will turn out to be good, in 




1 The charged lepton and Dirac-neutrino Yukawa couplings themselves receive small corrections induced by A° e ; 
the effect of going in the basis where charged leptons are diagonal is accounted for by a redefinition of the c's in 
EqM. 
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all the subsequent discussion we drop the O-superscript above the parameters of and U°. For 
later convenience, notice also that 
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(ii) 



4 General Implications of MLFV for thermal leptogenesis 

The MLFV scenario described above has six independent CP-violating phases. As shown in 
Appendix [B]they can be characterized in terms of weak-basis invariants that in turn are in direct 
correspondence with linearly independent combinations of the asymmetries e{ a . The analysis of 
weak-basis invariants (see Appendix [B] for more details) allows us to draw the following general 
conclusions on the viability of leptogenesis both in the one-flavour regime > 10 12 GeV) as 
well as in the flavour regime (M® < 10 12 GeV): 

• In the one-flavour regime, in order to produce the lepton asymmetry one needs non-zero 

(2) 

second-order RH mass splitting (i.e. c\ ^ 0) and H ^ I, i.e. CP violation (CPV) in the 
RH sector ED. 



• Flavour effects open at least in principle two new regimes for MLFV-leptogenesis which are 
not allowed in the one- flavour limit: (i) the case in which RH mass splitting is induced only 
(or mainly) by c^'. This situation requires H ^ I, namely CPV in the RH sector; (ii) the 
case in which CPV arises only from PMNS phases |23| . namely H = I and M v is real. In 

(2) 

this case the lepton asymmetry is proportional to c\ . 

We first study the Yukawa-induced RH mass splitting. The largest correction term is naturally 
expected to be the one associated to cW in Eq. (JM, so that many considerations about the RH 
neutrino spectrum and the frequency of the oscillations in m^'-*' can already be drawn. This 
program is straightforward in the limit of small (pi, which from the numerical analysis will turn 
out to be the relevant one in order to achieve viable leptogenesis. In this regime one can then 
expand the relevant matrices in series of <pi. In particular, one has from Eq.Q 

e~ 1 ^ \/m^ci2C 2 3 -e~ l ^ A /m2Ci 2 .S23 +£'{ ( t > i), (12) 



/m 2 Ci2C 2 3 -e 2 ^/m 2 Ci 2 S23 
/ rh~3S 2 3 V^3 C 23 

where the O-superscript is understood for rrii and the parameters of U°. In addition 

mi i(f>i(mi + m 2 ) ifaimi + iris) 

mi +771,2) "72 i(f>?j{m2 H 

-i(p2(mi + 7773) -7^3(7772 + mz) 7773 




M° 

K~-y%-\ -i<t>i{mi + m 2 ) 777 2 i4>^{m2+m^) | + 0(<pi(j)j) (13) 



and 

hp ~t~ h v — () x 

yZ 
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1m\ (p2<p3{ m i + m 2 + 2m 3 ) -<j)i^{m\ + 2m 2 + m 3 ) 

x| 0203(^1 + 1^2 + 2m 3 ) 2m 2 (/>i</> 2 (2mi + m 2 + m 3 ) 

-0i(/>3(mi + 2m 2 + m 3 ) 0i^ 2 (2mi + m 2 + m 3 ) 2m 3 

+ (14) 
The leading correction for RH neutrino masses is simply linked to the light neutrino mass spectrum 



M v = Ml 



1 + 2c (i)^^ (1 + CWi)) 



(15) 



Let us consider NH and IH for light neutrinos in turn. From the form of M v it turns out that in 
first approximation, 

/0 1\ 

U NH ss / , [7 /H « 1 J = J . (16) 

\0 1 0/ 

Aiming at a compact notation, we introduce the matrix m = diag(m;, rrtj, Wft,), where m; < < 
m/i (light, intermediate, heavy); clearly, for NH Z = l,i = 2, h = 3 while for IH/ = 3, i = l,h = 2. 



Then, defining Am, = rrij — rrij one has from Eq.(15): 

Mi » Ml (l + 2c« ^) , AM, » 2c«Am, (^Y , (17) 



So we have identified the parameters that control the frequency of the oscillating see Eq.Q. 

Moreover, the RH neutrino widths and the amplitude A„ of the sin-term in mA ij ' can then be 
expressed as 



m , 



r , w ( ^ ) , A j% = PS ^r-r^_ • (18) 

J 8vr V v J ' 3 AM^ 16vr cW Am, V ; 

We are now in a position to study the condition for the importance of quantum effects in the 
CP asymmetries, namely AMji < IV Using the expressions above, in the MLFV framework this 
condition can be cast into 

16tt c (1) J - < 1 , (19 

mi 

implying that the difference between the classical and quantum approach will increase as 
becomes smaller and/or the pair m^-mi becomes degenerate. In the next sections we will study 
numerically and analytically the dependence of the baryon asymmetry on cW and the lightest 
neutrino mass mi, controlling the degeneracy of the light neutrino spectrum. 

Both in the classical and quantum case, in order for e<J' to be non vanishing, it is crucial to 
have non diagonal entries in h u = = Uh^W . Since U is a small rotation, one has h v ~ h®. 



From Eq.(13) it turns out that the simplest situation in which h v is non-diagonal occurs when at 
least one among the <$>% is non vanishing, namely H ^ 1. On the other hand, in the case H = 1, 
small real non diagonal elements for h u have to be generated as an effect of U; in this case the CP 
violation needed for a non-zero asymmetry arises from the PMNS phases in {\ v )i a . These two 
situations are going to be discussed separately in the next sections. 

Before turning to analyse the differences between the quantum and classical approaches fo- 
cussing on these two relevant cases, we show the expressions for the washout parameters. For 
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both NH and IH, the washout parameters K a = £V Ki a and K\ = J2 a Kia (see appendix [A] for 
more details) are given by: 

K e = (micf 2 + m 2 s? 2 + m 3 s i3) I ' m * > (20) 

= (miC2 3 S? 2 + m 2C? 2 C23 + 7713523)/^* , ( 21 ) 
if T = (miSi 2 S2 3 + m 2 Ci2si3 + m 3C23)/ m * , ( 22 ) 
= rrii/m* (» = 1,2,3) . (23) 

It is clear that none of the K a {a = e, (j,, r) can be small: the weak-washout regime cannot 
be obtained in the framework of MLFV. However, let us emphasise that for NH the potentially 
smallest washout parameter is K e . On the contrary, for IH all the K a are comparable and 
relatively large (in particular, for small m 3 one has the relation K e ~ mi/m* ~ 2K fl ^ T ). 

In what follows, we shall focus our attention on the full flavoured regima^J M® < 10 9 GeV, 
and we shall illustrate the impact of quantum effects in the two cases in which CP asymmetries 
arise either from the RH {H 7^ I) sector or from the PMNS phases (H = I). 



5 CP Violation from H ^ I 



We first analyse the situation c^ 1 ** 3> and H 7^ /, in which the lepton asymmetries are propor- 
tional to the angles (pi. For NH, let us focus on the asymmetries associated to the least washed-out 
flavour, that is the e-flavour. Introducing the notation = l/(8ir)(gs + S« and by 

working at first order in 0j, we find: 



m 2 mi + m 2 (21) 

tie = \ = Ci 2 Si2<?>l-F^ ' -'cos 
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S12S13(P3 
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where fh = v 2 /M®. For zero CPV in PMNS the ej Q are maximised. Notice that for some values 
of the PMNS phases, one can suppress tjb] however, if all fa 7^ 0, there is no choice of PMNS 
phases that allows all three e to vanish^} 



(24) 



5.1 Dependences on c^ 1 ' and mi 

In order to explore numerically the dependence on we set cl = and choose M® = 10 9 
GeV. Then, we switch on in turn each cj>i = 10 -4 , setting to zero the two others. Because of the 
smallness of the fa, we will have 777^ ~ 777^ and U ~ U°; for definiteness, we select the parameters 
of U° to be: 9 23 = vr/4, 12 = 35°, 13 = 10~ 3 , 5 = a\ = a 2 = 0. 

For 777; = 10~ 3 eV (thick red line) and 777/ = 1CP 1 eV (thin blue line), the upper plots of Figs, [l] 
and [2] show the ratio of the asymmetry computed including quantum effects (t?^ u ) to the classical 

2 For M° > 10 9 GeV the results do not change significantly. 

3 Indeed, the term proportional to {4>i, <j>2, 4>'a\ vanishes when {0:2 — ai,2S + a±,2S + 02} — it + mod(27r), 
respectively. The latter two conditions imply Q2 — cti = + mod(27r), in conflict with the first condition. 
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Figure 1: Dependence on in normal hierarchy. Thick red and thin blue lines correspond 
to mi = 1CT 3 eV and mi = 10 1 eV, respectively. Top: The absolute value of the ratio of the 
baryon asymmetry with quantum effects (f?£ U ) and without quantum effects {rfg). Bottom: The 
absolute values of rf^ (solid lines) and rfg (dashed lines). 




Figure 2: Dependence on in inverse hierarchy. Thick red and thin blue lines correspond to 
mi = 10" 3 eV and mi = 10 1 eV, respectively. Top: The absolute value of the ratio of the 
baryon asymmetry with quantum effects (f?^ u ) and without quantum effects (r]g). Bottom: The 
absolute values of rfg (solid lines) and r/^ (dashed lines). 
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Figure 3: Dependence on mi, in normal hierarchy. Thick red and thin blue lines correspond to 
cW = 6 x 1CP 3 and = 2 x 1CP 5 , respectively. Top: The absolute value of the ratio of the 
baryon asymmetry with quantum effects (f?£ U ) and without quantum effects {rfg). Bottom: The 
absolute values of rf^ (solid lines) and rfg (dashed lines). 




-6 -5 -4 -3 -2 -1 -6 -5 -4 -3 -2 -1 -6 -5 -4 -3 -2 -1 

Log 10 [m 3 /eV] Log l0 [m 3 /eV] Log 10 [m_,/eV] 




Figure 4: Dependence on m.3, in inverse hierarchy. Thick red and thin blue lines correspond to 
cW = 6 x 10 -3 and = 2 x 10 -5 , respectively. Top: The absolute value of the ratio of the 
baryon asymmetry with quantum effects (f?£ U ) and without quantum effects (rig). Bottom: The 
absolute values of rfg (solid lines) and rf^ (dashed lines). 
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asymmetry (?7r), for NH and IH respectively. The corresponding lower plots show the absolute 
value of t]b with (solid line) and without (dashed line) quantum effects. In the classical case, 
the resonant point is usually reached for ~ 0.01, the only exception being IH with fa ^ 
where it is reached for = 0(1). Let us recall that the natural range for radiatively-induced 
c^ 1 ) is around l/(47r) 2 = 6 x 10 -3 . As can be seen, quantum effects induce an enhancement of t]b 
at small values of Such an enhancement is caused mainly by the second (sin) term of the 
function For both NH and IH quantum effects seem to be important for < 10~ 3 , the 

only exception being the case of nearly degenerate light neutrinos with fa as the main source of 
CPV in leptogenesis. In the latter case, values of up to 10 -2 could do the job. 

The figures also show that with degenerate light neutrinos quantum effects are larger, as 
expected, but t]b is more suppressed. Since the baryon asymmetry scales linearly with fa one 
would need these phases to be larger than about 10~ 4 to reproduce the observed value of the 
baryon asymmetry. The dependence on mi is shown in more detail by the plots of Figs. [3] and [4] 
for NH and IH respectively, selecting two representative values of c^>: = 6 x 10~ 3 (radiatively 
induced) and = 2 x 10~ 5 . Notice that the inclusion of the terms proportional to \ setting 
for instance = (c*- 1 )) 2 in Eq.([6|, does not change the above results. 



5.2 Analytical estimates 



In this subsection we will derive a simple analytical estimate for the ratio 



R 



m\ 



\cl\ 



(25) 



and check that it is in agreement with the numerical findings. 

The dominant contribution to the asymmetry will come from the least washed-out flavour /3, 
e.g. the one with the smallest Kp. Therefore, as a first approximation we take: 



R ~ 



\qu| 



(26) 



cM->0 



In general, Yaa receives three contributions from the out-of-equilibrium decays of the three iVj's. 
Each of these terms is proportional to the CP asymmetry e^g, which contains (X v \t)ijm^\ with 
j 7^ i. The expression (13) for X^Xt tells us that if e.g. only fa ^ 0, the only non-zero off-diagonal 



entries are those with = (1,2) and (2, 1). 

In the limit c^ 1 ' — > 0, the mass splittings AMji become small (compared to the scale M°), 
and so the arguments of the periodic functions in m^'i' are also small, in the range of z where 
the lepton asymmetry is generated. The sin 2 term in m™' is completely negligible; instead, the 
sin term has in front an amplitude proportional to 1/AMji. Therefore, the dependence on 
cancels in this limit and the resulting quantity is 



m 



(*) 



2H(Mx) 



(27) 



This disappearance of explains the plateaux observed in the upper plots of Figs. [I] and [2] in 
the limit of very small 
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Let us consider the case where only fa ^ 0. Looking at the expression (13), we only have 
two non-zero asymmetries, say ejp and ej-p, but we choose the one containing the least wash- 
out parameter Kj, since it corresponds to the most out-of-equilibrium RH neutrino whose decay 
produces the largest asymmetry; so we just take tjp{z) ~ —ZjpKjZ 2 /2, where ejp is the CP 
asymmetry without the memory factor. 

Our choices for the lightest neutrino masses (mi = 10~ 3 eV, 1CT 1 eV) make all of the K^s 
greater than or of order 1 . In such a case, it is a good approximation to take dY^ jdz ~ dY^/dz = 
— z 2 K,i(z) / 4g* (here and in the following K,\,K,2 are the modified Bessel functions of the first and 
second kind, respectively). Then, an approximate solution to the Boltzmann equations (Eqs. (61 ) 
and (62) in Appendix |A|) is found using the steepest descent method: 

o^r 1 - -±r rtee^KMe-^zr^w) (28) 

(29) 



4#* 

e J/3 1 



9* Kf 3 \A l 3 l3 \z 1 

where Z\ satisfies the condition: K ^^^ z\K,\(z{) + — jg^f^ — 0- A good analytical approxi- 
mation for zi, valid for Kp > 1, is z\ ~ 3.47 [log(Kp)] ^ 4 . 

Taking into account the quantum correction factor (|27l) gives: 



(YA,r - ~ T dze^K^e-^^^'^') (30 ) 



4^* 



8g* 



Kj I dzz 4 /Ci(z)e~^f^J z «*"Mi*j ( 31 ) 



2g* K/3\A/3{3 



z 2 • (32) 



where satisfies the condition: K ^^^ z\K\(z-i) + -J^ — ~ 0. Since one expects zi^ > 1, 

they satisfy almost the same condition; thus it is reasonable to take: Z\ ~ z,2 = The expressions 
(29) and (32) lead to the simple estimate: 

i?(0, ^ 0) ~ ^f 2 ~ 10 Kj [log(^l^l)] 1 - 28 , (33) 

where we neglected 0(1) factors, since we are only interested in getting the order of magnitude 
of the importance of the quantum effects over the classical approximation. If we had taken lower 
values of mi, thus giving Kj = min(Xj) < 1, the analytical estimate would have proceeded in a 
different way. Though, R receives a leading contribution still proportional to Kj as in (33), plus 
corrections going like Jj-^, where z eq is defined as the "time" at which the number density of 
Nj reaches the equilibrium one: Y^.(z eq ) = Y^(z eq ). For the reference values we considered, this 

additional term is ~ 10 -2 Kn/Kj. As described later, such a correction might be important in 
the limit of very small Kj, i.e. very small mi, since the leading term is suppressed in that limit. 

Let us now apply the analytical result (33) to the specific cases analysed numerically in the 
previous subsection. 



Normal Hierarchy. With the parameters chosen for the numerical analysis, the least washed- 
out flavour is /3 = e, while the Ki's are in the order K\ < K2 < K3. Thus, the previous arguments 
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lead to 



12(01 ^ 0) ~ 12(0 2 ± 0) ~ 10 #1 [logCfgAeel)] 1 - 28 , 12(0 3 / 0) ~ 10 K 2 [logCKel^l)] 1 - 28 . (34) 

In the case m\ = 10~ 3 eV, one has K\ ~ 0.9, K2 — 8.4 and !T e ~ 3.4, giving 

12(0i ^ 0) ~ 12(0 2 / 0) ~ O(10) , 12(0 3 / 0) ~ 0(1O 2 ) . (35) 

In the case m\ = 10 _1 eV, one has K\ ~ 93.5, K 2 ^ 93.8 and K e ~ 93.6, giving 

12(01 + 0) ~ 12(0 2 ^ 0) ~ 12(03 / 0) ~ O(10 4 ) . (36) 

Such estimates are in good agreement with the numerical values for rfg /rf^ in the c^ 1 ' — ► limit, 
shown in the upper plots of Fig. [TJ 

Inverse Hierarchy. In this case the least washed-out flavour is /3 = fj,, while the IQ's are in 
the order K3 < K\ < K 2 . Repeating the same analysis as before, one finds 

12(0! ^ 0) ~ 10 Id pogC^I^I)] 1 - 28 , 12(02 / 0) ~ 12(03 ^ 0) ~ 10 K 3 [log(^|A w |)] L28 . 

(37) 

In the case m® = 10 3 eV, one has K\ ~ 46, K3 ~ 0.9 and ~ 24, giving 

Rfa ^ 0) ~ O(10 3 ) , 12(02 + 0) ~ 12(03 ^ 0) ~ O(10 2 ) . (38) 

In the case 771° = 10 _1 eV, one has iTi ~ 104, K 2 ~ 93 and ~ 10 2 , giving 

12(0i + 0) ~ 12(02 / 0) ~ 12(03 / 0) ~ O(10 4 ) . (39) 

Again, these orders of magnitudes agree very well with the cW — ► limit found by numerical 
integration, as one can see from the upper plots of Fig. |2j 



As anticipated before, the same analytical approximation (33) qualitatively explains also part 
of the dependence on mj, for very small (the thin blue line in the upper plots of Figs. [3] and 
Indeed, by considering the various K a and Ki as functions of the light neutrino masses, as given 



by Eqs. (20)-(23), one can recover the decreasing behaviour of 12 as mi becomes small. In some 



cases, though, the main term in 12 given by (33) turns out to be proportional to mi itself. They 



are the cases 0i,02 7^ for NH and 02,03 7^ for IH. In such situations, when mi is too small 
(compared to m* ~ 10 -3 eV) the leading contribution to 12 is suppressed and the next-to-leading 
term proportional to 1/Kj starts being important. This describes the change of behaviour which 
can be seen in the upper plots of Figs. [3] and [4] in the cases we mentioned above. 

The fact that such a change in 12 occurs at different values of mi for NH and IH can be seen as 
follows. For NH, the next-to-leading order term contains K e /K\, as functions of the light neutrino 
masses; since S13 is tiny, the main contribution will come from the term with m 2 /mi = m^/ mi. On 
the contrary, for IH the largest term in the ratio K^/K^ is the one with m\/m 3 = m^/mi, which 
is generally an order of magnitude greater than mi/ 'mi for NH. Therefore, the next-to-leading 
term starts being relevant, in IH, at a larger value of mi than in NH. 

In this Section we have considered the effect of the splitting induced by - namely 0, as 
sources of CPV. We have found that what is important for quantum effects to be sizable is not the 
magnitude of the washout parameters K, which is independent of in the MLFV framework. 
Rather, what is crucial is the frequency and the amplitude of the sin oscillating term, which have 
to be respectively small and large. This requires small and/or strong degeneracy among the 
light neutrinos rrij and nij for the quantum effects to be sizeable. In the next Section we will 
analyse the case in which H = I. 
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6 CP violation from light neutrino mixing (H = I). 

MLFV-leptogenesis is viable in the limit H = I only in the "flavoured" regime, provided that 
the heavy neutrino mass matrix receives a radiative splitting proportional to the charged-lepton 
Yukawa couplings {5M V = c± M®(h e + /i<T)). Here we derive the analytic dependence of the 
asymmetries €i a on the parameters of the model. 

The CP asymmetries read 

Im (Aj,) iQ 

e* a = E pM — " . (4°) 

where, following the notation of Section 3, \ u = U\® and UM U U T = M v . In the limit U — ► 1, 
where A -> A° (see Eq. @), it is easy to verify that the €i a vanish. The key ingredient in obtaining 
a non-zero result is to have non-diagonal entries in h u = \ v \l and correspondingly a non trivial 
U matrix. In order to estimate the off-diagonal entries of h v and U we perform a perturbative 
expansion of M v assuming <C c^: 

M u = M° j/ + 2c« h° v + (2ci 2) + + 4 2) ) {hlf + cf (h e + hi) + • • • | , (41) 
= M + A . (42) 

Here M is a real diagonal matrix, with eigenvalues Afj — M^[l + 2c^nii/fh + (2c^ 2) + 4 2) + 
4 2 ^)(mj/m) 2 ], and A is a real symmetric matrix given by 

' c { ^M^{h e + h T e ) NH , 

A = <! (43) 
cf ) M^I(h e + hJ)I T IH, 



with the I defined in Eq. (16). M v is diagonalised perturbatively by the real orthogonal matrix 
U = I + T (T T = -T), with 

Tij = Mi — Mj ■ (44) 
To first order in c± /c^ the explicit epression for \„\t = Uh° u U^ is then 

A fm -m l r (2) f Re (^ NH ' 



m(M 4 - Mj) c« 



Re(Ih e I T )ij IH . 



Using the above result in the expression of 6j a (along with (A°)j Q (A°)^), one obtains for NH: 

A2) 



,7Vff 

rW 

while for IH: 



2) „ r ^ 

d=e,a,T j^i K ) 



4? = -C E (^) 2 E E ^ W« fe^) Im (C/* m C/ Qn ) 1 . (47) 

C „ \ %J / I //(/ I 



i> / * — ' * — ' i m 

P=e,fi,T j^i m,n 
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In the above expressions the dominant term arises when j3 = r in the sum over charged lepton 
flavours. Keeping only f3 = r and defining 



the asymmetries read: 



NH 



NH 
e 2a 



-3a 



.(2) 



V v J I rh 12 m 13 

+ CT fTUry r F (l ) 2)^ $ (a)_ jF (3,2)^ $ (a) 

C W V -u / L m 12 m 16 



c(i) 

.(2) 



+ - 



.(2) 
-A 

:(l) 



F (l,3)!^*(«) +i r(2,3)!^ $ (a) 

m d m ^ 



(48) 

(49) 
(50) 
(51) 



The IH case is obtained by a straightforward permutation of indices in the m,- and 3>f.- factors. 



We are now in a position to identify the dependence of the baryon asymmetry on the low-energy 
CP violating phases, contained in the factors . In the case of NH, only are relevant, as 
the e flavour is the least washed-out one: 



(e) _ 
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S 12 C 12 S 23 C 13 sin («l - «2) + 0(4?) 



(e) 
13 

(e) 
23 



S12C12S23C23C? 3 S13 COS ( sin ^ + ^ + 0(s? 3 ), 
S12C12S23C23C? 3 S13 COS sin + 5^ + 0(s 



Is) 



(52) 
(53) 
(54) 



For IH, all flavours are roughly equally washed out, so one needs also &[j' T ^ ■ Expanding to first 
non-trivial order in si 3 we find: 



(m) 

12 — 
M _ 



1 



2222 ■ / 

s 12Cl2S23 c 23 sm («l 



02) 



13 

(M) 
23 

(r) 
12 

(r) 
13 

to 
23 



2 2 2 2 
s 12 s 23 c 23 c 13 smQ; l 



C?2«23 c 23 c 13 sin «2 



s i2 c i2 s 23 sin(ai - a 2 ) , 



2 2 2 2 
s 12«23 c 23 c 13 sm «l 



C?2 S 23 C 23 C 13 sin «2 



(55) 
(56) 
(57) 
(58) 
(59) 
(60) 



The above expressions explain well the two main qualitative features of this scenario: 



• In order to have a lepton asymmetry which does not vanish in the limit si 3 
necessary to have non-vanishing Majorana phases. 



it is 



• In the limit a.\ = a 2 = the lepton asymmetry is proportional to S13 sin 8 but the overall 
scale is substantially smaller than in the generic case with non-vanishing Majorana phases. 
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Figure 5: Left: tjb vs. («i — 02) /2 varying mi < 10 2 eV, S13 < 0.1, and < 5 < 180° in the 
NH scenario. Right: t]b vs. 5 for ct\ = 02 = 0, setting mi = 10 -3 eV and S13 = 0.03 in the NH 
scenario. In both plots = 10~ 2 , c^ 1 ) = 10 _1 , M„ = 10 9 GeV, and red (blue) points correspond 
to the inclusion (exclusion) of quantum effects. 




Figure 6: Dependence on mi, in normal hierarchy and H = 1. Thick red and thin blue lines 
correspond to cW = 6 x 10~ 3 and = 2 x 10 -5 , respectively; the other RH mass splitting 
coefficients are taken to be: c[ \ 3 = an d = 0.1c*- 1 ). The PMNS phases and S13 are chosen to 
be ot\ = 36°, «2 = 60°, 5 = 0, s%3 = 10~ 3 and a± = 02 = 0, 5 = 7r/4, S13 = 0.1. Top; The absolute 
value of the ratio of the baryon asymmetry with quantum effects (rj B u ) and without quantum 
effects (t/^). Bottom: The absolute values of rj B u (solid lines) and t] b (dashed lines). 

Qualitatively, these two observations agree with the conclusions obtained in Ref. |23j by means 
of a numerical study0 They are also qualitatively confirmed by our numerical evaluation of the 

4 A precise numerical comparison with Ref. [23] is not possible given the different assumptions about the structure 
of the cf. 
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Figure 7: Dependence on 7723, in inverse hierarchy. Thick red and thin blue lines correspond 
to = 6 x 10~ 3 and = 2 x 10 -5 , respectively; the other RH mass splitting coefficients 
are taken to be: c{\ 3 = and cjj 2 = O.lc^. The PMNS phases and S13 are chosen to be 
a.\ = 36°, a2 = 60°, 6 = 0, S13 = and a.\ = «2 = 0,(5 = 7r/4,si3 = 0.1. Top: The absolute 
value of the ratio of the baryon asymmetry with quantum effects (^ u ) and without quantum 
effects (t/g). Bottom: The absolute values of rf^ (solid lines) and rfg (dashed lines). 

baryon asymmetry taking into account quantum effects, as shown in Fig. [5] 

It should be stressed that the above observations are only the main qualitative features and 
do not describe in detail all the allowed possibilities. In particular, the complete numerical study 
shows that the rjs oc si3sin<5 behaviour in the ot\ = ai = limit, derived from a linear expansion 
in S13 and in the approximation K e <C if» Tl is a good approximation only for S13 <C 0.1 and 
small mi. For large S13 and large m\ values the 0(sf 3 ) terms generate a sizable correction to 
the formally leading si3sin5 term. As a result, the asymmetry assumes a form of the type 
r)B oc sin 5(1 + fecosJ), with b = 0(1). 

As far as the size of the baryon asymmetry is concerned, if the Majorana phases are large 
one can easily reach values of r\s which are consistent with the experimental observations (see 
Fig.|5]left). This is not possible if the Majorana phases are set to zero, at least in the single-Higgs 
scenario considered so far. In this case t\b turns out to be smaller than the experimental value by 
at least one order of magnitude. However, the overall normalization of t]b changes if we consider 
a two-Higgs doublet scenario (such as the Higgs sector of the MSSM), where the charged-lepton 
Yukawa couplings are enhanced by tan/? = (Hjj) /(Hp) > 1. In this case t\b oc tan 2 /3 and the 
overall normalization can easily reach the experimental value even if the Majorana phases are set 
to zero. The MLFV framework with H = I, vanishing Majorana phases and large tan /3, provides 
a concrete scenario where the observed baryon asymmetry is directly linked to the measurable 
CP violating phase of the low-energy neutrino mass matrix. 

We conclude that the impact of quantum effects in the H = I case is qualitatively very similar 
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to the H 7^ I scenario: memory effects become sizable for small splitting (in particular small 
c^) and in the limit of light degenerate neutrinos. An illustration of the impact of these effects 
is provided by the the plots in Fig. |6j|7j As can be noted, the effects can be quite dramatic, 
especially for IH. However, it should also be stressed that the largest relative impact is obtained 
in regions of the parameter space where t]b is well below the experimental value. 



7 Conclusions 



In this paper we have studied resonant leptogenesis in the MLFV framework where it is assumed 
that the charged-lepton and the neutrino Yukawa couplings are the only irreducible sources of 
lepton-flavour symmetry breaking. In such a framework, the heavy RH neutrinos are highly de- 
generate in mass and their decays in the early Universe may give rise to the observed baryon 
asymmetry through the mechanism of resonant leptogenesis. Previous studies on the viability 
of leptogenesis in the MLFV framework |21| |22| [23] have been based on the assumption that 
classical Boltzmann equations suffice to analyse the dynamical generation of the baryon asymme- 
try. However, it has been recently shown |24[ [26] that quantum Boltzmann equations are a more 
appropriate tool to study such a dynamics when the heavy RH neutrinos are degenerate in mass. 
Indeed, the quantum Boltzmann equations obtained starting from the non-equilibrium quantum 
field theory reveal that the CP asymmetry is a time-dependent oscillatory function which reduces 
to the value obtained in the classical approach only if the oscillation time is much larger than the 
interaction time. In resonant leptogensis this is not the case. 

We have shown both analytically and numerically that neglecting the time dependence of the 
CP asymmetry may underestimate the baryon asymmetry by several orders of magnitude when 
a strong degeneracy among heavy RH neutrinos and small mass splittings in the light neutrino 
sectors are present. This is true both when the CP phases come from the RH sector (phases in 
the matrix H) and when they come entirely from the left-handed sector {H = I) and may be 
identified with the low energy PMNS phases. 
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A Conventions for Boltzmann equations 

In the early Universe the quantum numbers conserved by sphaleron interactions are the A a = 
B/3 — L a . The pair of Boltzmann equations describing the generation of the baryon asymmetry 
are 



dY Ni 
dz 



-Di(Y Ni -Y%), (61) 
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= - Y, eiaDi " Y n) - W « I A ™ I Y Aa , (62) 

i 

where at equilibrium the N{ number density normalised to the entropy density of the universe is 
Y^ q = zf K^te) I (4<?*) , where Zi = z^/xl, g* = 106.75 and X^te) is a modified Bessel function of 
the second kind. The washout parameters are generically defined by 

_ T(N t -> £ a H) 

H(T = Mi) ' lMJ 

and we also make use of the quantities: K a = Yli^ia and Ki = ^2 a Ki a . In terms of the 
parameters of the model, Ki a = \(\ u )i a \ 2 v 2 / (Mim*) with to* m 10~ 3 eV. The decay and washout 
terms appearing in the Boltzmann equations are 

A = Ki x t z , W a = V ~K iQ v 7 ^ Kite) ^ 3 , (64) 

K-2teJ t-r 1 4 



while the matrix A is given bjjj 

A = -diag(151/179, 344/537, 344/537), for M x < 10 9 GeV. (65) 

For 10 9 GeV < Mi < 10 12 GeV and M x > 10 12 GeV, the two-flavour and the one-flavour regimes 
should be applied, respectively [2"9"] . 

Finally, 

a 

to be compared with the measured value Yb = (8.7 ± 0.3) x 10 -11 or with the baryon asymmetry 
normalized with respect to the photon number density (instead of the entropy density) tjb = 
(6.3 ±0.3) x 10~ 10 . 



B Analysis of CP violating weak-basis invariants 



The independent CP-violating phases of the model can be characterised in terms of weak-basis 
invariants, i.e. quantities that are insensitive to changes of basis or re-phasing of the lepton fields. 
The MLFV scenario under investigation has six independent CPV invariants coming from the 
Yukawa sector. 

The simplest necessary conditions for CP invariance can be cast in the following weak-basis 
invariant form 



Bi 

B 2 
B 3 



ImTr 
ImTr 
ImTr 



hu {MlM v )M* h* v M, 
h v (MtM u ) 2 M* h* v M v 







o , 



h v (MtM u ) 2 M* h* v M v {MlM v 



o , 



(67) 
(68) 
(69) 



3 We used the approximation in which A is a diagonal matrix. For the full expression see refs. |27l 1281 125] . 
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where h v = A v At and M v denotes a generic heavy neutrino mass term. The invariants -8^2,3 are 
independent and survive in the limit A e —* 0. One can construct three other independent invariants 
that explicitly [^ involve A e by simply replacing one h u entry in Bi^.3 with h e = X u X\X e Xl: 

Bi = ImTr \h„ (M}M V )MZ h*M v ] = , (70) 
B 2 = ImTr \h v (MlM u ) 2 M* 7i*Mj,J = , (71) 
fc„ {MlM v ) 2 M* v H* e M v {MlM v 



B, 



ImTr 







(72) 



Note that Bi 2 3 and Bi 2 3 are in direct correspondence with (linearly independent combi- 
nations of) the CP asymmetries e« relevant for flavoured leptogenesis. In particular, Bi.,2,3 
correspond to the combinations relevant for the 1-flavour regime. This can be seen by working in 
the weak-basis where M„ and A e are diagonal. In this basis, for example, B\ and B\ read: 



i?! = £ MiM^Mf - M?) Im [(A^ a (At ) Qi (A.At 

i<j a=e,iM,T 

2 

Bi = J] MiMjiM] '- Mf) ^ Im [(A,) ia (At ) QJ (A.At ) 



2<j 



(73) 



If M v is proportional to the identity, the hermiticity of h v>e and the cyclic property of the 
trace operation imply that the Sj, Bi vanish identically. The next step is to break the degeneracy 
of the heavy neutrinos in a way consistent with the MLFV hypothesis. Selecting for convenience 
the 0- superscript basis, let us now investigate under which conditions on M v of Eq.Q and \® of 
Eq.([8| the Bj,Bj are non- vanishing. 

• If we confine ourselves to terms quadratic in the Yukawa couplings, i.e. proportional to 
then we have that: 

— Bi = because of the hermiticity of h® and properties of the trace operator. This 
implies that in the unflavoured regime the CP asymmetries vanish; 

— Bi / as long as H ^ I. Since Bi oc c^ImTr [(/i°) a /i*] = (a is some integer) and 
the restriction H = I implies that h„ is real, it turns out that leptogenesis is possible 
in the flavoured regime only if H 7^ I. 

• If we include those terms in Mr that are quartic in the Yukawa couplings, i.e. proportional 

(2) 

to c\ 2 3 4, then: 

— Bi 7^ if any of the c*p 7^ 0, as long as H / /. If H = I then not only h® is real, 
but also M v and therefore Bx2,3 = 0. So one concludes that with H = I unflavoured 
leptogenesis is not viable; 

(2) 

— Bi ^ if any of the c\ 7^ 0. In this case, even setting H = I leads to non-zero Bi as 

(2) 

long as c\ 7= 0. So one concludes that leptogenesis with H = / is potentially viable 

(2) 

only in the flavoured regime with c\ ' + 0. 



Of course A e can also appear implicitly through the contributions to M v allowed by MLFV. 
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In conclusion, flavour effects open at least in principle two new regimes for MLFV-leptogenesis 
which are not allowed in the 1-flavour limit: (i) the case in which RH mass splitting is induced 
only (or mainly) by This situation requires H ^ I, namely CPV in the RH sector; (ii) the 
case in which CPV arises only from PMNS phases, namely H = I and M v is real. 
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